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Cosmological backreaction suggests a link between structure formation and the expansion history 
of the Universe. In order to quantitatively examine this connection, we dynamically investigate 
a volume partition of the Universe into over- and underdense regions. This allows us to trace 
structure formation using the volume fraction of the overdense regions Xm as its characterizing 
parameter. Employing results from cosmological perturbation theory and extrapolating the leading 
mode into the nonlinear regime, we construct a three-parameter model for the effective cosmic 
expansion history, involving Xmq, the matter density and the Hubble rate H-Dq of today's 

Universe. Taking standard values for and H-d^ as well as a reasonable value for Xmq, that 
we derive from A^-body simulations, we determine the corresponding amounts of backreaction and 
spatial curvature. We find that the obtained values that are sufficient to generate today's structure 
also lead to a ACDM-like behavior of the scale factor, parametrized by the same parameters 
and Ht>q, but without a cosmological constant. However, the temporal behavior of Xm does not 
faithfully reproduce the structure formation history. Surprisingly, however, the model matches with 
structure formation with the assumption of a low matter content, ~ 3%, a result that hints to 
a different interpretation of part of the backreaction effect as kinematical Dark Matter. 

A complementary investigation assumes the ACDM fit-model for the evolution of the global 
scale factor by imposing a global replacement of the cosmological constant through backreaction, 
and also supposes that a Newtonian simulation of structure formation provides the correct volume 
partition into over- and underdense regions. From these assumptions we derive the corresponding 
evolution laws for backreaction and spatial curvature on the partitioned domains. We find the 
correct scaling limit predicted by perturbation theory, which allows us to rederive higher-order 
results from perturbation theory on the evolution laws for backreaction and curvature analytically. 
This strong backreaction scenario can explain structure formation and Dark Energy simultaneously. 

We conclude that these results represent a conceptually appeaUng approach towards a solution of 
the Dark Energy and coincidence problems. Open problems are the still too large amplitude of initial 
perturbations that are required for the scenarios proposed, and the role of Dark Matter that may 
be partially taken by backreaction effects. Both drawbacks point to the need of a reinterpretation 
of observational data in the new framework. 

PACS numbers: 98.80.-k, 98.65.Dx, 95.35.+d, 95.36.+X, 98.80.Es, 98.80.Jk 



I. INTRODUCTION 

One decade has passed since the first long distance SN 
measurements (see [l|-[3| for the latest data) led to the be- 
lief in an accelerating expansion of our Universe and the 
postulation of Dark Energy (even if the effect is model 
dependent d, [H and there are some possibilities that do 
not need acceleration (H, 0|)- There is still no convinc- 
ing theoretical understanding of the experimental result, 
and the creativity of theoretical physicists has led to a 
plethora of proposals on the nature of this new energy 
component, see e.g. [8]. 

In the same decade, however, initially independent of 
the development in the Dark Energy sector, our knowl- 
edge about another, long-standing problem in relativistic 
cosmology has grown to an extent that allowed to connect 
it to the Dark Energy problem: the question of how ef- 



fectively inhomogeneities in matter and geometry would 
influence the global behavior of the Universe [9] . As gen- 
eral relativity is far too complex to simply solve the Ein- 
stein equations for an inhomogeneous matter distribution 
and calculate the effect through tensorial averaging tech- 
niques, there have been several approaches to get a han- 
dle on this problem |10|. The one which lies at the basis 
of this paper was introduced in and uses a truncated 
hierarchy of the averaged scalar parts of Einstein's equa- 
tions. As will be recalled in Section III A( this strategy 
reduces the problem to a set of three coupled equations 
for four independent variables characterizing the state of 
some spatial domain in the Universe, its volume, mean 
density, mean scalar curvature and backreaction, linked 
to the degree of inhomogeneity of the domain's deforma- 
tion. For details see Sec. Ill Al and fTll-[T3|. 

Based on this framework the connection of the two 



2 



above-mentioned problems was made in 116| and [15|, 
where the authors proposed that the backreaction of the 
inhomogeneities developing in the era of structure for- 
mation could lead to the effect interpreted as accelerated 
expansion. This would also solve the coincidence prob- 
lem, i.e. answer the question why the acceleration takes 
place at a moment of the evolution when the Universe 
leaves its homogeneous initial state. For this and other 
more theoretically founded reasons, there have already 
been some applications of the formalism of the averaged 
equations to structure formation. One has been carried 
out by Rasanen |17] (generalizing his earlier work on a 
two-scale model [18]). He used a model of structure for- 
mation in which the extreme points of the initial over- 
density field are evolved separately using a spherical LTB 
model. His treatment did not give evidence for an accel- 
erated expansion, but he could identify a signal of onset 
of an eventually strong backreaction regime at the right 
time. Another study that also follows a multiscale ap- 
proach as in the present paper has been performed by 
Wiltshire [l^, [lOl • He separated the universe model into 
"voids" and "walls" and studied the backreaction effects 
resulting from fluctuations due to this partitioning. To 
solve the equations he assumed the backreaction term to 
vanish on each of the subregions and set the curvature on 
the overdense regions to zero (as in the two-scale model 
of Rasanen [18]). For the under dense regions he consid- 
ered a Friedmann-like constant curvature term. Under 
these assumptions he also did not find volume accelera- 
tion, but described an effect that should result from the 
different lapse of time in underdense void regions and 
overdense matter dominated regions. This allowed him 
to fit his model to a fair number of data [l9|, [U, [22| , 
despite his assumption that the backreaction effect on 
the individual domains was absent. We shall come back 
to his approach later. Other interesting considerations 
about a multiscale approach to self-gravitating systems 
may be found in [23]. 

We wish to take one step back and return to the more 
general problem (i.e. without using Rasanen's or Wilt- 
shire's assumptions), by presenting the general frame- 
work for separations of the averaged equations into sub- 
domains of the Universe. This has been begun in j23| by 
investigating a general volume partition of cosmological 
hyper surfaces. We will show that this separation is con- 
sistently possible also for the evolution equations. The 
advantage of not neglecting any effect is impaired by the 
caveat not to be able to close the system of averaged 
equations without further assumptions on evolution laws 
for backreaction and curvature. We are going to inves- 
tigate a volume partition of the Universe into initially 
overdense A4-regions and underdense ^-regions in anal- 
ogy with Wiltshire's approach. However, instead of set- 
ting the backreaction on A4 and £ to zero, we will first 
use perturbative results on the early evolution of back- 
reaction [26i-28]: in perturbation theory one finds 
that the backreaction term on a subdomain of a matter 
dominated universe model decays only as with the 



scale factor (a result found in various papers in relativis- 
tic cosmology [l5|, [IB, l26l - [28| that confirms the formally 
analogous situation for the backreaction term in Newto- 
nian cosmology [25^). The latter work also shows that it 
is the leading term in a nonperturbative approximation 
based on the exact averaged equations and the Zel'dovich 
approximation p9l-[32| for the fluctuations. Furthermore, 
this work selects the leading term as the dominant con- 
tribution on large, "quasilinear" scales. Note here that 
in some work this leading mode is dismissed due to the 
property that its coefficient function is a full divergence 
and must hence vanish by imposing periodic boundary 
conditions. While this latter remark is true the 
mode is nevertheless there in the interior of a periodic 
Newtonian or quasi-Newtonian simulation box. More- 
over, in a proper relativist ic theory these terms are not 
full divergences and there are no boundary conditions 
to be satisfied apart from the constraints on the initial 
hyper surface. 

As a result of this behavior of the leading mode the 
importance of the backreaction term, although being a 
second-order quantity, quickly rises with respect to the 
matter density which goes as a~^. This perturbative 
mode falls on an exact scaling solution (3^ of the gen- 
eral problem and we, hence, consider evolution laws for 
backreaction and curvature by extrapolation of the per- 
turbative mode along this particular scaling solution. We 
emphasize that this extrapolation of the leading pertur- 
bative mode already furnishes a nonperturbative model 
due to the fact that the mode is referred to an evolv- 
ing background that includes backreaction (sH, [36| , (g^I • 
As an aside we here note that a perturbation theory on 
the evolving background given by the averaged equations 
has not yet been concisely investigated; forthcoming work 
will especially focus on the formulation of such a theory, 
which will be key to an adequate generalization of per- 
turbation theories in the averaging framework. 

Our second approach uses the general formalism and 
is specialized to a concrete model by making use of N- 
body simulation data and the assumption that backre- 
action globally acts like a cosmological constant. This 
latter approach illustrates a strong backreaction scenario 
that simultaneously describes accelerated expansion and 
the structure formation history. TV-body simulation data 
have recently also been used to estimate the backreaction 
in a quasi-Newtonian setting through an estimation of 
the post-Newtonian potential [38]. We here allow for a 
more general reinterpretation of Newtonian simulations 
that - combined with the relativistic equations - takes 
scalar curvature into account. 

This paper is organized as follows: Section [III shortly 
reviews the averaged equations and the occurrence of the 
backreaction terms. Then the relevant separation for- 
mulae are presented and the consistency of the resulting 
evolution equations is checked. The question of how ac- 
celerated expansion may be understood in the context of 
these equations is addressed in Subsection [Aj Section ITTTl 
uses the assumption of the a~ ^-scaling of the backreac- 
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tion terms to clarify the connection between structures 
today, the initial density fluctuations and the acceler- 
ated expansion. The comparison of llllDl shows that the 
specific model is not sufficient to explain structure for- 
mation and points out why this could also not have been 
expected. The discussion will concentrate on how one has 
to devise the evolution of inhomogeneities in the context 
of the averaged equations in order to fully explain the 
phenomenon of Dark Energy. 

Section IIVI investigates how the perturbative a~^- 
behavior has to be extended to combine structure for- 
mation and accelerated expansion. In IVK we motivate 
the choice of the a~ ^-scaling for the backreaction term 
on M and 8 and rederive the results of [25l|. f27l. [2l| in 
our context. We also perform a quantitative estimation 
of the amount of initial backreaction necessary for struc- 
ture formation. Finally, we give quantitative conclusions 
in IVB. and discuss in Sec. IVII what extensions could be 
made to the models considered and comment on the need 
for a substantial reinterpretation of observational data. 

II. VOLUME PARTITIONS OF THE AVERAGE 
UNIVERSE 

A. The averaged Einstein equations for compact 
spatial domains 

The basis for our investigation of the link of structure 
formation to the expansion history of the Universe is to 
consider integral properties of scalar variables in a pres- 
sureless - i.e. dust - universe model, given in [11], later 
refined to be applicable to perfect fiuid matter models 
in [l^, where also compact formulations of the averaged 
equations of [ll| can be found. The space time is fo- 
liated into fiow-orthogonal hypersurfaces featuring the 
line-element 

ds^ = -df + QijdX'dX^ , (1) 

where the proper time t labels the hypersurfaces and 
are Gaussian normal coordinates (locating free-falling 
fiuid elements or generalized fundamental observers) in 
the hypersurfaces. gij is the full inhomogeneous three- 
metric of the hypersurfaces of constant proper time. Such 
a split corresponds in the ADM formalism to the choices 
of a constant lapse and a vanishing shift. Extensions of 
this formalism to nonvanishing shift and tilted foliations 
have been given in (40| and (39|, in (4l| with applications 
to perturbation theory. For considerations of covariance 
and gauge issues see [42]. 

On these hypersurfaces we want to study the evolution 
of compact spatial domains P, comoving with the fiuid. 
This latter property ensures that the domain is frozen 
into the general three-metric, i.e. its shape encodes the 
geometrical structure of the local inhomogeneities. Note 
that in Newtonian spacetimes it can be shown that the 
morphology of the boundary of a comoving domain al- 
ready contains comprehensive information of the mat- 



ter distribution, including all higher-order correlations. 
The morphology-characterizing set of integral geometri- 
cal measures is known as the Minkowski functionals, from 
which the backreaction term can be built (see [l3| for de- 
tails). One fundamental quantity characterizing such a 
domain is its volume, which is the only such measure 
used here, since we wish to address questions related to 
the size of the domains and its time-derivatives only, 

■■= ^ dn, , (2) 

where d/ig := ^ (3)g{t, X^X"^, X^)dX^ dX^ dX^ . From 
the domain's volume one may define a scale factor 

..W=(p^j , (3) 

encoding the average stretch of all directions of the do- 
main. For wild changes of the shape of the initial domain 
Vi one might want to know more about the evolution of 
other morphological characteristics to deduce directional 
expansion information, and would therefore have to ex- 
tend the analysis. 

Concentrating on the volume and the effective scale 
factor alone, one can derive, from the Einstein equations 
with a pressureless ffuid source, the following set of equa- 
tions governing its evolution: 

3^ = -47rG(^?)p + Qp+A (4) 

ml = 8nG (g)^ - \ {7^)^ - ^Qp + A (5) 
= dt {q)j, + 3ffp {q)^ , (6) 
where the average over scalar quantities is defined as 

and where ^, IZ and H^) denote the local matter density, 
the Ricci scalar of the three-metric gij^ and the domain 
dependent Hubble rate Ht> '= dv/dv^ respectively. The 
kinematical backreaction is defined as 

2x>:=^((^%-W^)-2(a^)^, (8) 

where 6 is the local expansion rate and := l/2(Jij(j'^^ 
is the squared rate of shear. Note that Hj) is defined as 
Ht> = 1/3 (6) j^. Qt> is composed of the variance of the 
local expansion rates, (^^)^ — {^)\)^ and the averaged 
shear scalar (cr^)^^ on the domain under consideration. 
For a homogeneous domain it is zero. It therefore encodes 
the departure from homogeneity and is supposed to be 
particularly important in the late, inhomogeneous phase 
of the Universe at the epoch of structure formation. 

The integrability condition connecting Eqs. dU and (0 
reads 

a^^dt {al{n)v) = -a^'S, {a%Qv) , (9) 
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which already shows an important feature of the averaged 
equations as they in general couple the evolution of the 
backreaction term, and hence extrinsic curvature inho- 
mogeneities (or in this picture matter inhomogeneities), 
to the average intrinsic curvature. Unlike in the case of 
a standard Friedmannian evolution the curvature term 
here is not restricted to an scaling behavior but is 
dynamical in the sense that it may be an arbitrary func- 
tion of ax)' We emphasize that the essential effect of 
backreaction models is not a large magnitude of Qx>^ but 
a dynamical coupling of a nonvanishing Qx) to the aver- 
aged scalar curvature, changing the temporal behavior of 
this latter. 

One may also define cosmological parameters on the 
domain investigated. In analogy to the Friedmannian 
case they are derived from the Hamilton constraint (|5]) by 
a division by ZR^- different spatial domains, which 
will be denoted by P, JV[ and we formulate the cos- 
mological parameters in a generic way by taking T to 
denote one of the domains P, A^, f (this will hold for 
the whole text, but depending on the context T may also 
only denote M. and £). The definitions are 



SttG 



(0) 



in) 



T 



ml 



A 

ml 



(10) 



where the decision to divide by Rj^ for every T will be- 
come clear from the definition of A4, f and V in the next 
section. Using those definitions, the Hamilton constraint 
([5]) for a domain T reads 



^2 



(11) 



which means that the dimensionless parameters only 
add up to 1 for the domain V. On other domains they 
may add up to a more or less arbitrary positive value, de- 
pending on whether the corresponding region F expands 
faster or slower than the V region. 

To point out the analogy with the /c = Friedmann 
equations one may recast (|1])-(|6|) combining the backreac- 
tion and the curvature terms to an effective density and 
pressure: 

3^ = A - 47rG(pfff + 3p?ff) ; 31f|, = A + SttG^^Sj 
= ^ffF + 3ffp(f?fff+pSr) , (12) 
where the effective densities are defined as 



1 



1 



V 

Peff • = 



{^)v (13) 



WttG 



In this sense, Qx) and (71)^ may be combined to some 
kind of dark fluid component that is commonly referred 



to as X-matter. One quantity characterizing this X- 
matter is its equation of state 



^A,efr 



^eff 



ie) 



V 



(14) 



which is an effective one due to the fact that backreaction 
and curvature give rise to an effective energy density and 
pressure. 

B. Separation formulas for arbitrary partitions 

After the short review of the averaging formalism we 
present in the following how the resulting equations can 
be separated, if one wants to consider the behavior of 
subdomains of some "global" region P, which we may 
assume to be associated with a (postulated) scale of ho- 
mogeneity. We consider, a priori, arbitrary partitions of 
the spatial hypersurfaces into subregions J^^, which them- 
selves consist of elementary space entities T^'^^ that may 
be associated with some averaging length scale. The idea 
will be to choose entities ^j"^ in a way that they share 
some common properties. In our case this property will 
be the value of the overdensity field, but the separation 
is general and also valid for other choices. To be more 
precise we want to have V = U^Ti where := U^J-'j"^ 

and T^^""^ n = for all a ^ f3 and £ ^ m. In the se- 
quel we follow closely the investigation in [24] of a general 
volume partitioning of cosmological hypersurfaces. 

The average of the scalar valued function / on the 
domain V, 



if) 



V 



\V\ 



fdUg 



(15) 



defined in (O, where \V\g denotes its volume \V\g := 
fj^dfig, may then be split into the averages of / on the 
subregions J^^, in the form 



(16) 



where we have introduced 



\v\ 



(17) 



i.e. the volume fraction of the subregion T^. Regard- 
ing the dynamical equations we wish to separate, i.e. dU, 
([5]) and ([9]), Equation (fT6|) directly provides the expres- 
sions for the scalar functions ^, IZ and Hx '•= 1/3 {0)jy. 
Only Qx as defined in ([8j) does not split in this simple 
way due to the (0)'^- 
provides 



term. A short calculation rather 



= ^ A,Q, + 3 ^ A, Am {He - H^f 



(18) 
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as the correct formula. Qi is defined as in (|8|) with in- 
stead of V. The shear part {o''^)jr is completely absorbed 
in Q^, whereas the variance of the local expansion rates, 
(6>^)^ — (^)^, is partly contained in Qi but also gener- 
ates the extra term 3^^_^^ A^A^ {H£ — Hm)^ - This is 
because the part of the variance that is present in Q^, 
namely (^^)jr^ — ^^^^ takes into account points 

inside F^. To restore the variance that comes from com- 
bining points of Fi with others in J^rn-, the extra term 
containing the averaged Hubble rates emerges. 

One may also define a scale factor ai analogous to (|3|) 
for each of the subregions J^^. As their definition implies 
that they are disjoint, it follows that = \^^\g 
and we may therefore define = a|. An important 
detail when using the equation in this form is that, at the 
initial time when ax) is equal to one, the scale factors 
of the subregions will not be equal to one as well. This 
different normalization has to be taken into account when 
scaling the local M and £ parameters. 

Differentiating a^y = ^^a^ twice with respect to the 
foliation time and using the result for d^, we finally pro- 
vide the relation that links the acceleration of the scale 
factors of the subdomains to the global one: 



- ^ A^A^ {H£ — Hmf 



(19) 



As an immediate consequence one can see that even when 
the subregions decelerate, the second term of ([T9|) may 
counterbalance the first one to lead to global accelerated 
expansion. We will examine this property, which appears 
generically in averaged models due to correlations of local 
expansion rates, in Appendix [Al 

In the following we will be mainly considering the case 
where we have only two types of subregions Tf. over- and 
under dense regions. We will define them at some initial 
time and call M. T\ the part of the hypersurfaces that 
consist of all elementary regions A^(") := jr(") with an 
initial overdensity and 8 := the complementary part 
£ := P n A4, i.e. the one with the initially underdense 
regions. The formulae ([T6)) and (p!8]) then simplify to give 



Hx) = Hm + (1 — ^m) Hg 



(20) 



with the same expression vahd also for (g)^ and {TZ)x>j 
and 



Qv = \mQm + Ci- - >^m) Qs 

+QXM{l-^M){HM-H£f 



(21) 



In both cases we used A^ = 1, which results from the 
fact that the are disjoint, and defined \m •= 1-^1 / 1^1 
and A^: := |f | / This general separation formula al- 
lows us to illustrate the simplifications applied by Wilt- 
shire in his two-scale model |19| . As he assumes the over- 
dense Al-regions to have no curvature, by Eq. ([9|) the 
only possible nonzero Q;\4-term would have a strongly 
decaying a^-behavior and may therefore equally well be 



set to zero. For the underdense f-regions the situation 
is similar as he assumes them to have the Friedmann- 
like a^^ constant curvature term. This again makes the 
left-hand side of Eq. (|9]) vanish and leads to the choice 
of Qs = 0. This reduces the general formula (|2T]) to its 
third term on the right-hand side, which corresponds to 
his Eq. (31) in fl^^ that he derived for the case of vanish- 
ing shear. This is consistent, as we remarked below ([T8|) . 
because the shear term of Qx) is completely contained 
in Qm and Qg. Considering the overall formalism his 
choice is conservative in the sense that one would indeed 
expect the biggest backreaction effect to stem from the 
variance of the expansion rates between the A4- and £- 
regions and not from those on A4 and £ themselves, as 
they are chosen to have similar properties and therefore 
similar expansion rates. On the other hand, setting the 
local backreaction to zero restricts the new interesting 
feature of the averaging formalism, i.e. the coupling of 
backreaction to curvature, to the P-regions only. Fur- 
thermore, as we argue in Subsection IIII C\ it may be el- 
igible to allow for shear on the -regions, which may 
(and will) imply a negative Qm (i-e., a Dark Matter be- 
havior over J\4). We will therefore not restrict ourselves 
to the assumptions of constant curvature, but use other 
well-motivated conditions in Section ITTTl 



C. Consistent split of the dynamical equations 

With the knowledge of the relation between the quanti- 
ties on the subregions with those on the global domain P, 
one may now ask, how the separation affects the evolu- 
tion equations for those quantities. Therefore, we insert 
the expressions ([T6|) for Hx^ {q)x {^)x 
Qx into dl]), ([5]) and (|9]). A straightforward calculation 
shows that the equations take the following form: 



= 

= 
= 



£ 



Qe + 



-AttG {g)^ + + A - 3 



ai (t) 



(22) 
(23) 



Y,>^4^f9t{aj{n),)^aj'dt{a',Qi)] , (24) 



i.e. they may be split into a sum in which the equations 
on the subregions take the same form as those of the 
global region (H])-® and their contribution is weighted 
with the volume fraction of the respective subregion. As 
Equations (|4])-([6]) hold for an arbitrary domain V and 
therefore also for the subregions J^^, Equations (|22|) - ([24|) 
show that the separation advocated in the previous sec- 
tion is also consistent at the level of the evolution equa- 
tions. This is not surprising as the separation procedure 
is straightforward and the equations are supposed to hold 
on any domain, but was, especially in view of the nonlin- 
ear form of (p!8|) . not entirely clear when just looking at 
the formulas. The consistent split assures that, if we have 
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found a solution for the quantities on the subregions Ti 
and use the relations of the previous section to calculate 
those on the global domain, then we will automatically 
obtain a solution of the averaged Equations (|4])-(|6]) on 
this global domain. 



III. BACKREACTION SCENARIO BASED ON 
EXTRAPOLATING THE LEADING 
PERTURBATIVE MODE 

After setting up the equations that link different sub- 
regions of a certain partitioning of space to the whole in 
the preceding section, we will now have a closer look at 
the simplest partitioning. As already mentioned at the 
end of IIIB( we are interested in a subdivision into over- 
and under dense regions. These will be labeled by M. and 
f , respectively, and the classification is made at some ini- 
tial time and then kept fixed. As the averaged equations 
are derived for a certain region Tn^ that is assumed to 
keep its identity, an exchange of elementary space enti- 
ties T^^^ between the two classes of regions M. and E 
would spoil the applicability of (H])-® on M. and ^, and 
therefore would lead to some complicated modified form 
which would introduce even more unknown parameters. 
This is the reason for keeping the identification of Tg- 
regions, even if this will introduce some problems with 
the clear cut interpretation of M. and for initially over- 
dense regions may become underdense in the course of 
evolution. A refined description that helps to evade this 
problem will be mentioned in the discussion in Sec. IVII 

The motivation of distinguishing between the hA- and 
f-classes of regions is to get an explicit handle on struc- 
ture formation. Implicitly it is already present in the av- 
eraging formalism mainly via the {(O'^)^ — (6>)^^-term 

in the kinematical backreaction defined by (|8j). As struc- 
ture formation implies a different evolution of the local 
expansion rate at different places in the Universe, this 
variance is expected to grow. The distinction now of- 
fers the possibility to examine the development of parts 
that are in the focus of experimental interest, namely 
"voids" (see e.g. [43] and references therein) and clusters 
of galaxies. Moreover, it allows us to trace the process of 
structure formation by the simplest parameter character- 
izing its history, i.e. the volume fraction of the overdense 
regions. The current section presents the outcome of this 
investigation for the special case of exact scaling solu- 
tions. 



A. Exact scaling solutions of the averaged 
equations 

As we have seen in Section III A| the averaged equa- 
tions do not form a closed system. One closure strategy 
is to impose a specific equation of state. This is simi- 
lar to closing the Friedmannian equations, but here the 



equation of state is dynamical. As done in quintessence 
models for the X-matter content inserted as source into 
the standard Friedmann equations, the simplest assump- 
tion is to choose a constant equation of state parameter 
w. This means that we have a single scaling law for the 
corresponding matter component. If we proceed in a sim- 
ilar way for the backreaction and curvature terms in the 
averaged equations, we may choose (following [34]): 

Qt = Q^.a^jr ; (7^)^ = n^y^ , (25) 

where T stands for or ^ and where Ti indicates the 
initial state of the domain J^. By the integrability condi- 
tion ([9]), Qjr and (7^)jr are related. Inserting the above 
ansatz leads to two types of solutions: first, there is a 
solution with n = — 6 and p = —2. This one is not 
very interesting as this is the only case where the generic 
coupling of backreaction and curvature is absent. This 
case would correspond to fluctuation histories that are 
decoupled from a constant-curvature (Friedmann-like) 
behavior. Second, ([9]) will be satisfied if n = p. This is 
the case we are going to study in the following. It implies 
that Qjr and {1Z)-p are related by a constant r*^ 

Q^ = r^(7^)^ = r^7^^,aJ, (26) 
which is determined by (|9]) to be 



The effective equation of state for the X-matter (backre- 
action and curvature) component, defined in is for 
scaling solutions simply 

<eff = -^(« + 3). (28) 

For the explicit solution of the equations on M and E 
in this paragraph we will further specify n to —1. The 
motivation for this choice is, that Li and Schwarz [26| 
found, in a calculation using second-order perturbation 
theory, that the backreaction term may be expressed as 
a Laurent series with a leading term going like a^^. We 
will rederive this behavior, under similar conditions as 
in their case, in Section |Vl There we also argue that 
the leading term of the scalar curvature, going like 
is expected to vanish for plausible initial conditions. We 
will show that extrapolating the -behavior up to to- 
day runs into trouble and that higher terms in the series 
are expected to make an important contribution at the 
epoch of structure formation. Nevertheless, the a^^ scal- 
ing solution provides a first and simple example of why 
one would expect that backreaction through structure 
formation might cause accelerated expansion [44]. 

We may argue that we did not have to rely on assump- 
tions on the effective equation of state, if we exploited 
the knowledge of structure formation further. Such an 
attempt runs into problems as we do not have an inho- 
mogeneous relativistic metric model at our disposal that 
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would cover the details of nonlinear structure formation. 
On large scales we probably understand what is happen- 
ing, but we are already using this information by the 
known scaling behavior on these scales. In the late 
Universe with density contrasts of order 1, the closest 
to a realistic model at present is the one by Rasanen 
But also this has to rely on some unrealistic as- 
sumptions. We therefore try to derive the equation of 
state from the N-body simulation of structure formation 
in our second approach in Section IIV( where we use the 
knowledge on nonlinear structure formation encoded in 
the volume fraction of the simulation. 



B. Free parameters and constraints 

In order to derive the evolution of the global domain 
V under the assumption of the -scaling behavior for 
M. and we have to solve Equation ([5]) (the other inde- 
pendent Equation ([9|) is already satisfied by the scaling 
ansatz) which simplifies to 



' 4 



,^0 

7^Q 



for out of M and ^, where we have defined ^^^q := 
{^'^ + Iiistead of using and as^ as free pa- 

rameters we may equally well work with A^o ^Vq- 
To specify these and the other free parameters, namely 
^X)o' ^Mo^ ciVo and ^^^Q - with the index denot- 
ing the value of today - we have to make some further 
assumptions. One is that the density fluctuations in the 
Early Universe around z = 1000, where we choose our ini- 
tial time ti^ were Gaussian and very small as indicated by 
the results of the CMB experiments. For our parameters 
this means that A^j. ^ 0.5 as Gaussian fluctuations im- 
ply that there were as ma ny ov er- as underdense regions 
and we set a^. = a^;. = y^1/2. The assumption of being 
close to homogeneity implies {Q)vi ^ {Q)Mi ^ {Q)si' 1^ 
then follows that 



^0 (^Vo 



3H; 



-Do ^To 



(30) 



■J-Q 



which allows us to simplify Equation (|29]) by replacing 
Vt-^"" and Vt^ with . To reduce the set of unknown 
parameters even further, we may use (|29]) today which 
gives 



Mo 
1ZQ 



^0 



(31) 

and may be used together with (|2Q|) to eliminate ^^q- 
Another consistency condition finally allows us to get rid 
of the parameter as well, enabling us to fix our 

model without having to know anything about the values 



of the backreaction or curvature terms. The argument is 
as follows: consider the integral of Equation 



-Do 



a 



/1/2 



(32) 

which gives two functions tj\/[ (a>() and is, (ag) depending 
on the parameters in the integral. This may be stated as 

tj^iaj^) = tjr {aj^-.Vt-^Q.Hvo,^^^ ,avo,\M^' As there 
is only one time even in the separated equations, which 
is the time of the foliation, Im and ts have to be equal. 
(Note that here we do not introduce a phenomenological 
lapse between times of different regions as done by Wilt- 
shire, since this would require another (very involved) im- 
plementation of a multiscale foliation; we so must imply 
that Wiltshire's effect is negligible, but we do not claim 
that this is indeed the case.) Using this condition today, 
i.e. demanding that Im (<^A4o) — {(^So) gives another 
relation linking and to the rest of the param- 
eters. We may therefore eliminate as well, which 
leaves us with the four free parameters Hd^ , l^^o , a-Do 
and Xmo • ciVo and Hx)q only determine the scaling of the 
axis so that the parameters responsible for the shape of 
the curve of time evolution of aM {t) and ag {t) are 1]^° 
and Xmo- 



C. A specific numerical example 

In order to get some quantitative statements we will 
specify the remaining parameters to some plausible val- 
ues. For this is the one indicated by the experiments 
that established the cosmological concordance model. 
Even if it is not a priori clear that the domain depen- 
dent of the averaged equations should have the same 
value as the Qrn parameter whose value is determined 
from the data under the assumption of a homogeneous 
and isotropic universe model, this is the only thing we 
can do before the analysis of observations in the context 
of the averaged model, begun in [45|, is completed. In 
this spirit we will use for A^q some reasonable value de- 
rived from a A/'-body simulation of large-scale structure 
formation. These are in our context not the best sources 
of data either, as they are Newtonian and therefore do 
not include scalar curvature being a major player in the 
averaged equations. The reason why we think that the 
value of the parameter Xmq is nevertheless derivable from 
them, is first, that they are designed in a way to match 
best today's Universe and secondly, that, for the calcula- 
tion of the volume, even a substantial curvature deviation 
from a Euclidean space does not change the value very 
much. This is because the calculation of volumes is per- 
formed with the metric and its spatial derivatives that 
are linked to curvature do not enter in the comoving fo- 
liation used here (see (46| for clarifying remarks on that 
point). 
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Figure 1: Plots of some quantities derived from the single scaling model with Xmq = 0-09 and = 0.27. Top left: Comparison 
of the evolution of the different scale factors of our model and the one in the ACDM-case. On the right the corresponding 
Hubble rates. The one on Ai goes to zero in the end whereas Ht> approaches Hg asymptotically. Ha and Ht) are very close to 
each other. The middle row shows the Q-parameters on V and S. Those for M are left at the bottom. Finally, on the right 
the behavior of the effective equation of state as defined in fT^ . The a^^-scaling implies that w^j^ff and w^^eS both have a 
value of — |. 



We will explain further details of our analysis of the N- 
body data in Appendix IB II The value we will be using 
in this section has been obtained by dividing the analysis 
volume into blocks with sidelength 5h~^Mpc and evalu- 
ating whether they are over- or underdense. This leads to 



Xmo ^ 0-09. Together with fl'^^ ^ 0.27 and a^^ ^ 1000 
we may numerically evaluate Equation ([29|) which results 
in the curves of Fig. [TJ The value of ^ 1000 is used 
to identify "today" in our plots. It can be interpreted 
in line with the standard model, if the global volume is 
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Figure 2: Left: Comparison of the luminosity distances of our models (based on the template metric of [45|). with the one of a 
flat ACDM model with h = 0.7 and Qm = 0.27. On top the model of Section HVl where we force the volume scale factor ax) to 
follow the ACDM evolution. Despite this assumption, the changing curvature affects the luminosity distance. The luminosity 
distances in our models show a significant feature at a redshift of around 1, when compared with the best fit ACDM model, 
which may be looked for in the SN data. The curve below is the model of the current section. For comparison we also included 
the luminosity distance of the best fit model of [i^. Because of a different Hubble rate of /i = 0.7854 it lies below the others 
from the beginning. This model does not significantly show the distinct feature of the other two models around a redshift of 1, 
due to the assumption of a single-scale cosmology. 

Right: Values of Clarkson's C-function [72] for the best fit model of [45| (top), the model of Section [Tvl (middle), and the 
model of this section (bottom). Recall that for every Friedmann model C{z) vanishes exactly on all scales and for all redshifts. 
For the inhomogeneous models shown in the plot, this function has a minimum which may serve as observational evidence for 
the effective cosmologies, as proposed in [45]. As our multiscale model shows, it is not even necessary to measure derivatives 
of distance, since the feature is already present in the distance. 



comparable, i.e. if perturbations of the metric remain 
small. In any case, substantial changes are expected for 
the spatial derivatives of the metric, and consequently for 
the time-derivatives of the scale factor. 



1. Evolution of the model universe 

The top left graph of Fig. [1] shows the evolution of the 
different scale factors with cosmic time. For comparison 
it also contains a plot of the scale factor aA {t) for a flat 
standard FLRW model with a matter content of Qrn = 
0.27 and 73% Dark Energy. One can see that its evolution 
is very similar to the one of the volume scale factor a-z) (t) 
of the global domain V throughout the whole history. 

That this apparent match is also expected to hold true 
for observables is suggested by the analysis in Fig. [2l 
where we compare the luminosity distances at a given 
redshift of the average model and the standard ACDM 
model. To calculate the distances we used the effective 
template metric of [45], which provides a prescription 
of how to incorporate the changes in the average scalar 
curvature for the calculation of distances. We see that 
for the region where we have ample SN data, i.e. around 
a redshift of 1, the difference between our model and the 
ACDM model is less than 5%. 

Considering the evolution of and ag we find that 
they behave as one would expect from their nature. The 
overdense regions M begin their evolution as the un- 



derdense ^-regions due to their assumed similar matter 
content. Then their density increases by gravitational in- 
stability and the expansion slows down. The underdense 
regions which contain most of the voids continue their 
expansion. This may also be seen from the Hubble rates 
in the plot on the right. The expansion rate of the over- 
dense regions drops down to zero rather rapidly whereas 
the one of the underdense regions slows down much less. 
The shrinking percentage of the overdense regions causes 
Hjy to bend up again implying acceleration. The differ- 
ence to the FLRW Hubble rate is not yet visible, but as 
the evolution proceeds the cosmological constant of the 
FLRW model will force to rise in the future whereas 
H-r, will join the evolution of Hg. 

2. Energy content 

The plot of the ^l-parameters in Fig. [1] is consistent 
with the conclusion drawn from Eq. ^ that has 
to be less than — to give accelerated expansion. In 
addition, the magnitude of Vt^ underlines the statement 
that curvature is an important player in our model. The 
combination of these two facts has an interesting con- 
sequence: if, even for a small but nonzero backreaction 
parameter, the average scalar curvature parameter is of 
order 1, and therefore the actual matter parameter Vl^ 
is small, the model may explain Dark Energy without a 
large amount of backreaction. A small larger than 
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— does, however, not explain an accelerated expan- 
sion, but this latter property may actually not be needed 
to explain observations. This remark is particularly in- 
teresting in view of the observation in Subsection IIII D 3l 

That curvature is in fact more important than sug- 
gested by the restricted constant curvature FLRW mod- 
els has been shown recently in j46|. The plots in the 
following row clarify its origin. It stems from the under- 
dense regions which are supposed to contain most of the 
large voids that have developed negative curvature in the 
course of evolution. The value for the matter parameter 
of about 0.15 makes clear that one may not think of £ 
being composed of voids only. From its origin as the un- 
derdense half of the initial universe volume, it should be 
clear that it also has to contain regions that have devel- 
oped over densities. 

For the overdense regions the situation is different in 
that in the course of their evolution, positive curvature 
emerges. This is what one would expect due to the fact 
that they contain many gravitationally bound systems 
like clusters. Another observation that coincides with 
the intuitive picture we have, is the fact, that the back- 
reaction contribution is negative [note again the sign in 
the definition (p!Q|) ]. According to definition (|8|), nega- 
tive backreaction means that the region over which we 
average is shear dominated. In view of the filamentary 
structure of the matter distribution of the Universe this 
is also what we would expect. 

When interpreting the l]-parameters of A4 and £ one 
should have in mind that they do not have to add up to 
unity but, due to the definition of the ft^s using Ht> and 
not Hm and Hg, wih add up to H^^/H^ and Hg/H^. 
This means for A4 that the sum at the end will be close 
to 0, whereas the value on £ approaches 1. 

Finally, the plot in the last row on the right shows 
the evolution of the effective equation of state for the 
X-Matter component. For Ai and £ the value is, ac- 
cording to ([28|) . simply —2/3. For the global domain V 
we have an evolution from this limiting value to —0.95 at 
the maximum of structure formation which then relaxes 
back to ^ 0.9 today. It is interesting to note that the 
value stays relatively constant throughout the period of 
structure formation. 

This observation, i.e. that the a^^-scaling model on 
the partitioned domains naturally leads to an approxi- 
mate cosmological constant behavior on the homogeneity 
scale, may also be seen from Fig. [3l There it is shown, 
how the backreaction and curvature terms deviate from 
their initial a^^-scaling when structure formation sets 
in. They stay approximately constant for quite a long 
period of the evolution before the ^-regions finally dom- 
inate and the backreaction and curvature terms begin to 
decay faster again. 

If it were easily possible to interpret the quantities of 
the averaged model in terms of standard model ones, one 
might well imagine the value of the equation of state 
today of Fig. [T]to lie within the 1 — a boundary derived 
by the WMAP team for a time-dependent dark energy 
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Figure 3: Plot of the evolution of Qt) and (71)^^ in terms of the 
global scale factor a-v- For comparison we added a line with 
a simple -scaling and one that is constant. Qt) and (71)^ 
are normalized by — 6if|). so that the values at the initial 
time represent Qq and Q^'. We appreciate that the back- 
reaction terms feature an approximate cosmological constant 
behavior on the homogeneity scale despite the assumption of 
-scaling on the partitioned domains. Physically, this re- 
sult can be attributed to the expansion variance between the 
subdomains and, hence, this latter is identified as the key 
effect to produce a global Dark Energy-like behavior of the 
backreaction terms. 

equation of state (i^l • We do not show the corresponding 
plot however, because the comparison is not expected to 
be very meaningful as the analysis uses the FLRW model 
as a prior. We again have to ask the reader to await a 
more reliable interpretation of the data in the context of 
the averaged model. 



3. Dark Matter 

The comparison of the l]-evolution of the A4- and £- 
regions in Fig. [1] illustrates nicely the property of aver- 
aged models that, depending on the domain under con- 
sideration, the effective energy content may vary. Seen 
in terms of X-matter this means that this component 
behaves in a more Dark Energy-like way when we con- 
sider only the f -regions, since the backreaction is pos- 
itive and drives acceleration. For the Al-regions, how- 
ever, the backreaction term is negative and slows down 
the evolution in the same manner as a Dark Matter com- 
ponent would act (iil. If backreaction is an important 
contribution in today's Universe, this might imply that 
parts of the effect known as Dark Matter could be due 
to a combination of the backreaction and curvature term 
on A4-regions. This would perhaps open the possibility 
to reinterpret Dark Matter models in which the particles 
decay and the Dark Matter responsible for the acoustic 
oscillations might be replaced through its cosmological 
effect by a rising X-matter component. 

It is to be expected that only parts of the effect at- 
tributed to Dark Matter is modeled by the effect of a 
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negative backreaction, since there are many independent 
hnes of evidence for the existence of Dark Matter. While 
part of the Dark Matter problem related to the expan- 
sion history might well be attributed to backreaction, 
this would only alter the Late Universe constraints. The 
relative abundances of Dark Matter to baryonic matter 
in the Early Universe, determined by the CMB, are not 
expected to change if we assume a near-homogeneous 
initial state of the Universe. All CMB constraints that 
stem from a combination of primordial information and 
constraints derived from observations of the Late Uni- 
verse, however, will be subject to changes in course of 
a reinterpretation of observational data. To show what 
conditions the CMB data really impose, and that will 
have to be satisfied also by a backreaction model, we re- 
fer to a recent paper |49i] that analyses the CMB in a 
model independent way. 

If we choose for the X-matter component the descrip- 
tion of [s^l in terms of an effective scalar field, the "mor- 
phon", we would have two effective potentials depend- 
ing on whether we are in an over- or underdense region. 
They are of the form 

um{sm) = -^Hl,n^{-j^jhm-'{sM) (33) 
us {ss) = -^HlQi (7^^) * sinh-^ {ss) , (34) 

where := ^^/^^, sjr (t) := VS^^jr (t) and 

ujr (s) := SttGUj^ {^). These potentials are simply the 
specialization of the general result (37) of [s^l to the 
-scaling and the different negative sign for Qm- The 
sin~ ^-behavior reflects the fact that formally the Ai- 
regions recollapse. In reality the description as a comov- 
ing perfect fluid will, however, break down before so that 
only the rising branch is expected to be physical. The val- 
ues of the parameters on Ai (more precisely the fact that 
(7^)^ is positive and n = — 1, see [34] for details) imply, 
that we have an effective phantom field for the overdense 
regions. In spite of the usual interpretation in terms of 
Dark Energy, the different signs lead to a positive effec- 
tive pressure and therefore to a decelerating component, 
which hence rather acts in a matterlike way. The form 
of the potential for the parameters of the model of Fig. [1] 
has been plotted in (50| • The approach to describe the ef- 
fect of inhomogeneities by an effective scalar field allows 
us to connect the results obtained by the quintessence 
and scalar field community to the backreaction formal- 
ism and to reinterpret their fields and potentials in terms 
of physical quantities, i.e. the parameters of the averaged 
model. This is potentially interesting as sinh^-potentials 
have been shown to be able to behave like a Dark Energy 
component, e.g. in (5l| or [52| [where their (cosh —1)^ po- 
tential is just a sinh^ one by cosh {2x) — 1 = 2 sinh^ {x)]. 
For more detailed information to possible quintessence 
potentials, see (53| and references therein. 

Note that scalar fields have also been used to model 
Dark Matter halos [U and that other potentials for a 



scalar field have been advocated to implement the cosmo- 
logical evolution of such scalar field Dark Matter (55l , [56| . 
Consequently, it has been proposed to unify these Dark 
Matter and Dark Energy scalar fields into a single de- 
scription e.g. in [57I and Questions related to unifi- 
cation have been also addressed in a recent work on em- 
ploying the Chaplygin gas as effective equation of state 
in the present context [59]. Note also that the common 
strangeness of scalar field models, i.e. requiring parti- 
cles with masses of order of 10~^^eV is not present in 
the morphon picture as there the notion of a "particle" is 
only an effective one. 



D. Comparison with other calculations 

In order to judge the quality of the model presented in 
the previous section, we will now compare its predictions 
with other calculations in the literature. 



1. Magnitude of the initial backreaction 

We first consider the perturbative results of jsH, |4l|. 
There it has been shown that, at the recombination 
epoch, the magnitude of the effective energy contribution 
of backreaction and curvature is expected to be of the or- 
der of 2 X 10~^ on the scale of the horizon for a Hubble 
rate of h = 0.7 that has to be compared to a cos- 
mological constant contribution in the standard model of 
3 X 10~^ (as is demanded in the strong backreaction sce- 
nario below). The corresponding term in our model (on 
the scale of the simulation box), + becomes, for 
the special parameters of the previous section (the l/ajr 
scaling scenario and at an initial time where aj:, = 1), 
8 X 10~^ compared to the matter density parameter. This 
comparison illustrates, that even if the Early Universe is 
in a near to homogeneous state, where one can apply 
perturbation theory, the Late Universe can look strongly 
different by the effect of the leading order perturbative 
a^^-mode alone. When higher order -growing- terms 
of their Laurent series play an important role, or when 
nonperturbative terms enter due to a substantial injec- 
tion of backreaction at the stage of nonlinear structure 
formation, the situation will even become more different. 
As this comparison shows, the initial amplitude needed 
is still too large in the 1/ajr scaling model. However, 
our value was measured on the scale of structure, and 
may still drop by going to the horizon scale, which would 
reduce the discrepancy. 

2. Structure formation in the concordance model 

The second comparison we want to present is the one 
for the evolution of Xm • 

In the corresponding figure. Fig. |4l the line shows 
Xm (^) as calculated from our model, where we have 
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Figure 4: Top left: Comparison of the evolution of the parameter Xm calculated from the scale factors shown in Fig. [T] (straight 
line), with the values derived from the A^-body simulation (dots). The graph shows that for the A^-body simulation the 
evolution of Xm is faster at the beginning and slower at the end compared to the one for our model. 

Top right: Comparison of Xm as calculated from Eq. ()35p for ^ 0.03 (straight lines) with the data (dots) derived from the 

A^-body simulation for grid sizes of = 5, 10, 15, 20, 30, 50h~^Mpc (from bottom to top). The only input used to generate 

these curves was the last point of the A/'-body data and the value of . 

Bottom left: Same plot as top right but with = 5, 10, 15h"^Mpc only and = 0.018. 

Bottom right: Same plot as top right but with R = 20, 30, 50h"^Mpc only and = 0.035. 



chosen to present the results in terms of z to provide 
a more intuitive picture of the actual evolution of struc- 
ture formation. The points result from our analysis of 
the A^-body simulation [gO^ using a simple separation 
into blocks described in Appendix IB II The figure shows 
a clear discrepancy between the Newtonian A^-body sim- 
ulation and our averaged model. This is not surprising 
as we are approximating the whole Laurent series of Qjr 
and {TVjjr by their leading terms only. We will see in Sec- 
tion [Vl how the higher terms are connected to the shape 
of Xm (^)- Of course it may also be that a genuine rela- 
tivistic simulation might add corrections to the values of 
Xm 5 but for the moment we have to work with the models 
we have. Therefore, we will present, in Section IIV| the 
characteristic features of a model that describes struc- 
ture formation and accelerated expansion in the present 
framework of averaged models in a consistent manner. 



3. Structure formation in general a-p^ -scaling models 



We finally want to report on an observation that might 
become of interest in the upcoming reinterpretation of 
the observational data in the backreaction context. As 
described in Appendix IB 1) we conducted the analysis of 
the A/'-body data using blocks of various side lengths. In 
addition to the grid size of 5h~^Mpc, being at the basis of 
the result of the previous section, we also used spacings 
of 10 — 50h~^Mpc. The resulting values for Xm (^) are 
shown as dots in the graph on the right-hand side of 
Fig. m As we have seen on the left-hand side of this 
same figure, the -scaling model using the concordance 
value of 1]^° ^ 0.27 is not able to reproduce the structure 
formation history as observed in the A/'-body simulation. 
We may however try to derive the evolution of Xm (^) 
for different choices of The differential equation 



determining Xm (^) may be derived from the definition 
of the latter and reads: 

^ . _ Hm {Cty, /orN 

6 Am -Livyo^v^AM) 

where Hj) = Xm^m + (1 - ^m) Hg and Hm {clv, ^m), 
respectively. Hg (a-z)^ Xm) may be found by replacing 
aM civ^M and ag ax) {I — Xm) in Equation (|29|) . 
Using the relations between the different parameters dis- 
cussed in Subsection llllBI we may numerically solve (j35j) 
after having specified the two remaining parameters 
and Xmo' The result for l^^o ^ 0.03 is shown by the 
lines in the graph on the right-hand side of Fig. HJ There 
we have fixed, for each curve Xm {^i^m^ i ^Mo)i ^Mo 
to the end value of the result from the TV-body simu- 
lation. It is interesting to notice, that the common value 
of ~ 0.03 allows us to describe the evolution on di- 
verse scales quite well and that the shape of the function 
as defined by Equation (|35]) meets the actual form with- 
out having to use the theoretically unmotivated fitting 
ansatz that will be presented in Subsection IV B[ The fit 
may even be improved if one allows for a weak scale de- 
pendence of 1^^°, as for ^ 0.018 the small scales are 
fitted in an exquisite manner and for ^ 0.035 this is 
true for the large scale evolution. 

If this is just a coincidence or if this points to a deeper 
physical interpretation of the backreaction effect in terms 
of its capability to unify Dark Energy and Dark Matter in 
one effective fiuid cannot be finally decided from this in- 
vestigation. It is clear that the backreaction terms imply 
the emergence of both dark components, however, in an 
entangled way that does not allow us to assign Dark Mat- 
ter to over-dense and Dark Energy to underdense regions 
uniquely due to the changing variance between the two 
types of regions. We may speculate that, in the course 
of a reinterpretation of observations in the context of the 
averaged equations, the role of Dark Matter in struc- 
ture formation on large and intermediate scales might be 
attributed to average curvature and backreaction com- 
ponents that effectively yield the same distribution of 
visible matter as CDM-models (but see the constraining 
remarks of Sec. IIII C 3p . But it may also be that a low 
value of 1]^° ^0.03 turns out to be incompatible with 
other observations. These questions will be addressed 
in upcoming work when testing the a^^-scaling model 
against observations as done for a simple a^^-model in 
(45j (with an improved template metric). Apart from 
observational issues the model has to be compared to 
N-body simulations with baryonic and neutrino content 
only in order to quantitatively support this necessarily 
speculative aspect of our model. A naive interpretation 
of the low-l]^° model also runs into the same difficulties 
of high initial amplitudes of perturbations as in the old 
studies of baryonic universe models. Furthermore, since 
independent evidence for Dark Matter exists 0|, espe- 
cially on small scales and unlike the situation for Dark 
Energy, we cannot conclude from this work that Dark 
Matter may be fully identified with backreaction effects. 
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But, our study clearly indicates the need to exploit this 
aspect of the backreaction effect. 



E. Discussion 

It may be helpful to hold in for a moment and put 
into perspective what we have seen in this section. It 
has been shown that in a model that possesses the a^^- 
limiting behavior of the backreaction components and 
that starts from almost homogeneous initial conditions 
as in the standard model, the simple fact that there is a 
lot of structure today in terms of volume dominance of 
devoid regions (manifesting itself in our model through a 
value of Xm around 0.1) implies some A-like accelerated 
expansion of the volume scale factor ap. Even though 
this model is too basic as is indicated by the discrepancy 
seen in Fig. 21 it clearly illustrates the way in which one 
would expect backreaction to act. (Here we do not want 
to over-emphasize the result that a pure baryonic matter 
content would even do the job.) To avoid the impression 
that we just replaced the mysterious Dark Energy com- 
ponent by some other mysterious component, i.e. the 
backreaction term, let us emphasize the different physi- 
cal situation. In the general relativistic framework that 
lies at the basis of the averaged equations the emergence 
of structure is associated with a geometrical deforma- 
tion of the underlying spacetime. In the comoving time- 
synchronous slicing chosen to describe the dust universe 
of III A( there is no movement of matter particles with 
respect to the space. All inhomogeneities emerge from 
the distortions of space itself. This results in the emer- 
gence of intrinsic curvature of the spatial hypersurface, 
reflected in the three Ricci scalar 7^, as well as extrin- 
sic curvature of this same hypersurface. This extrinsic 
curvature is what makes up the backreaction as it has 
been shown that it may be deflned instead of kinemat- 
ical quantities as in Eq. ([8]), in terms of invariants of 
the extrinsic scalar curvature (see [11]). To gain some 
intuition it may be helpful to see Equation ^ not in 
an active sense stating that there is some fundamental 
component that forces the scale factor to accelerate, but 
rather in a passive sense that it traces the complex evo- 
lution resulting from the full Einstein equations with the 
only active component being the perfect fluid dust mat- 
ter. Then one may picture the evolution to be governed 
by gravitational instability that causes inhomogeneities 
to grow, which manifest itself in extrinsic curvature show- 
ing up in the backreaction term, but also - and this is 
the key issue - in intrinsic curvature of space in which 
structures emerge. How acceleration in this context may 
be understood will be discussed in Appendix [Al 

In this general-relativistic picture one may also see 
that, even without a cosmological constant, we have some 
kind of "Dark Energy" in our model universe that now 
has a clear physical interpretation: It is the curvature 
energy of the spatial hypersurface, communicated by an 
effective potential energy in the morphon fleld correspon- 
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dence of backreaction jsj. This new kind of "Dark En- 
ergy" therefore necessarily emerges when structure forms 
and it is this "Dark Energy" the title of this paper is re- 
ferring to. 

Having emphasized the role of intrinsic curvature one 
cannot argue that this picture should yield the same re- 
sults as the Newton-inspired picture of structure forming 
in space and matter moving around. If the averaged in- 
trinsic curvature evolves differently as compared with a 
constant curvature model - and this is a generic outcome 
of the fact that the kinematical backreaction term and 
the averaged scalar curvature are coupled - then curva- 
ture plays a substantial role for the structure formation 
history, but also for the interpretation of observational 
data. While metrical deviations from a flat space may be 
small, derivatives of the metric can be substantial, as has 
recently been shown by estimating its magnitude in (46| . 
In addition the existence of a strong average scalar curva- 
ture term today is not excluded by the data. A study on 
light propagation in statistically homogeneous universes 
[62| has shown, that the position of the first acoustic peak 
in the CMB spectrum is consistent with a non-negligible 
amount of curvature if the expansion history of the scale 
factor ax) follows a A CDM-e volution (see also [63|). As 
this condition is fulfilled in our model, the occurrence of 
= 1 today is not problematic, if observational data 
are interpreted in the backreaction context. 



IV. MODELING STRUCTURE FORMATION 
AND ACCELERATED EXPANSION: STRONG 
BACKREACTION SCENARIO 



structure formation in a model using the partitioning in- 
troduced in Section [III The philosophy is that we want 
to show for which evolution of the regional backreaction 
and curvature terms the phenomenological global ACDM 
model, that provides a good fit to a large number of data, 
may be understood to be the result of the evolution of 
the physical variables extrinsic and intrinsic curvature. 

A. Evolution of parameters 

In order to construct the nonperturbative model dis- 
cussed above, we first have to translate the assump- 
tions of ACDM evolution and TV-body structure forma- 
tion into quantities of our model. The ACDM behav- 
ior means 17^ + 1^^ = ^Iriedmann ^ ^j^-^j^ rcsults in 
Qv {t) = A and (7^)^ = — 3A. A is determined by re- 
quirmg + 1]^° ^ 0.7. This is possible, since the cos- 
mological constant is a particular exact scaling solution 
of the averaged equations. To find Xm {clv) we analyzed 
the A/'-body simulation [60] using this time a separation 
into subvolumes based on a Voronoi tesselation of the 
simulation volume. The resulting data points for five dif- 
ferent redshifts are shown in Fig. [5l They are fitted with 
a functional ansatz of the form 




and the best-fit parameters providing the curve of Fig. [5] 
were aM ^ 191 and 13m ^ 419. 



After having explored the partitioning approach with 
the help of a simple -scaling solution in the previous 
section, which provided the result that this leading per- 
turbative mode alone is able to account for a A-like ac- 
celerated expansion if we insert what we know about the 
matter content and the structures of today's Universe, 
it may be interesting to learn how the higher terms in 
the Laurent series of Qjr and {IZ) j, have to contribute to 
give rise to the evolution of structures as derived from the 
A/'-body simulation. In other words we are searching for 
the nonperturbative behavior of Qjr and (7^)jr that will 
match the simulated structures. One would like to get 
a handle on the expansion behavior from simply tracing 
this evolution with the help of \m (^) but as the system 
will still not close one has to find another condition to 
come to results not being based on assumptions on the 
scaling of the unknown backreaction and curvature terms 
themselves. Even if, in principle, all the quantities that 
figure in the averaged equations are measurable, this is 
very difficult in practice and we do not yet have a precise 
idea how they evolve. For the backreaction and curvature 
terms, Ref. |24] provides strategies on how they could be 
determined, but there are no quantitative results yet. 

We will therefore only show here that the averaged 
equations are able to describe accelerated expansion and 



1. Results for the scaling of and {TV} j: 

The resulting scale factors of the domains are shown 
in the upper right panel of Fig. [5l Comparing them to 
the evolution of the scale factors for the -scaling it 
becomes clear, that a faster structure formation in the 
A/'-body case, which was found in Fig. 21 requires aM to 
slow down much earlier. The size of the -regions is 
therefore nearly constant throughout a long period. This 
is again qualitatively as expected, because the overdense 
regions virialize and decouple from the overall expansion. 
as evolves similar to the a^^-case, but finally has to take 
over the accelerated expansion and therefore has a limit- 
ing behavior of a cosmological constant evolution. This 
limiting behavior might not be realistic as it implies a 
constant variance of expansion rates which is assumed to 
shrink as the fastest expanding regions will dominate the 
volume in the late-time limit. If there is an intermediate 
state providing this evolution is speculative. 

The lower panels of Fig. [5] show the evolution of 
Qm {clm) resp. {11) M (^m) (left) and Qs (ag) resp. 
{71) g (as) (right). They were multiplied by aM resp. ag 
to point out the -behavior at the beginning of the 
evolution. It is interesting that this limiting behavior, 
which is - as already mentioned - in accord with the 
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Figure 5: Two partition model for a ACDM behavior of the V scale factor ax) and the structure formation inferred from the 
A/'-body simulation. Left on top the fit of Xm to the data points of the Voronoi evaluation. To its right the dependence of 
the scale factors on the foliation time. Bottom left, functional dependence of backreaction and scalar curvature on A4 on the 
corresponding scale factor aM- Because of the double logarithmic scaling only the absolute value is shown. (7^)^ is positive, 
Qm starts negative and gets positive at around aM ~ HO. Bottom right, the same plot for and {71)^ as a function of ag. 
{IZ) ^ is negative, Qs positive at the beginning and in the end, with a negative period between 110 and as ^ 500. In both 

cases the terms Qjr and (7^)jr have been multiplied by ajr to emphasize the initial a^^-limit. 



perturbative result of Li and Schwarz (26|, [28|, arises nat- 
urally in the present setup. The reason why this is the 
case will be explained in Section |Vl It should be noted 
that at the beginning, the curvature on Ai is positive, 
the backreaction negative. For £ it is just the opposite. 
Regarding the further evolution it can be seen that the 
backreaction changes sign (on £ even twice) and that Qjr 
and {TV) j: shrink faster than a^^. The reason for this be- 
havior can be seen in the shape of Xm in Fig. [5l First we 
need a strongly different behavior of M- and f-regions 
to assure that Xm shrinks rapidly. Therefore Qm is neg- 
ative and acts like matter to slow down the expansion 
of A^, whereas Qg is positive to lead to a faster growth. 
Then the growth on 8 is slowed down by a negative Qg^ 
whereas Qm becomes positive to counterbalance the de- 
celeration of the matter component to lead to the nearly 
constant part of the evolution of the -scale factor of 
Fig. [5] upper right. The linear rise for 8 at the end en- 
codes the cosmological constant behavior, meaning that 



Qs (cLg) = {^)s (^^) ~ const., which is reflected by 
an increasing line in this figure as we are multiplying it 
by as. 

2. Derived quantities 

Figure[6]finally shows the functional form of the dimen- 
sionless parameters. The upper left graph simply shows 
the evolution imposed by the required cosmological con- 
stant behavior of T) going from the matter dominated 
era = 1 to today's value of = 0.27. The rest, 
which is dubbed Dark Energy in the standard concor- 
dance model is now represented by our X-matter and 
split into backreaction and curvature. The f-regions in 
the upper right graph develop rapidly "strong" negative 
curvature, which is however not due to a rise in the (JZ)g- 
term (as may be seen from Fig. [5]), but to the faster 
decrease of the (^)£;-contribution. The evolution on M. 
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Figure 6: Dimensionless parameters for the model deduced from the fit of Fig.O Top left the evolution on P, on the right the 
one on ^, left at the bottom the same for M. Bottom right the effective equation of state of the Xjr -component as defined in 
The impUcations are discussed in Section llV A 21 



shows that one has to be careful in the interpretation of 
the r^-parameters. Due to the approximately constant 
scale factor a^^, which results in a nearly constant (^)^, 
the division by Hx> which is decreasing makes ft^ in- 
crease rather strongly. If one interprets SH^^/ (SttG) as 
the critical density, this increase in reflects the emer- 
gence of a strong density contrast between the Al-regions 
and the averaged universe model. 

Another interesting plot is shown in Fig. [6] in the lower 
right panel. It presents the effective equation of state of 
the X-matter component of our model composed of Qjr 
and (71) j: as Xjr := Qjf + ('^)jr- '^Aeff ~ ~^ reflects the 
imposed cosmological constant behavior on P, whereas 
the initial values on £ and A^, t^^i^ff ~ '^a eff ~ ~2/3, 
are related to the a^^-limit shown in Fig. [5l ap- 
proaches — 1 and w-ji^^^ rapidly goes to zero. This means 
that it is scaling like a matter contribution with aj^ in 
the end, but as the sign of Qm is positive, it acts as 
"repulsive matter", counterbalancing the deceleration of 
the matter component to lead to a nearly constant qm- 
Again, this is expected physically. While our model does 



not describe the relevant small-scale effects like velocity 
dispersion and rotation, effects that stabilize the virial- 
ized regions, the N-body simulation does describe these 
effects. The evolution of the Q>[-term thus rather re- 
flects the behavior of a more general backreaction term 
as expressed in [l3|. 



B. Cosmic phase space 

Another possibility to characterize the evolution of 
Qjr (aj^) and (71) jr {aj^) is to trace their solutions in a 
"phase space" introduced in [34]. Its dimension is two, i.e. 
there is room for many homogeneous, almost-isotropic ef- 
fective states, while the homogeneous-isotropic solutions 
just form a line in this space. Every curve in this space 
represents a solution of the averaged equations, while 
straight lines represent the class of scaling solutions. It 
is shown schematically on the right-hand side of Fig. W\ 
The different sections are exhaustively discussed in (33 |. 
We will repeat what is necessary for our purpose in the 
following. The coordinates of this space are chosen to be 
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Figure 7: "Cosmic phase space" of the solutions of the averaged equations. On the right the schematic partition discussed in 
|33 |. On the left the curves for M- (dotted) and 8 (dashed) regions derived from the A^-body simulations. It should be noted 
that they are shown here in the same plot only for economic reasons and actually live in two different phase spaces. 



the matter parameter and the deceleration parameter 
g*^, defined as 



(37) 



This latter was chosen instead of alone in order to 
have an additional intuitive meaning. In this space, every 
straight line passing through the center, which represents 
the EdS model, corresponds to an elementary scaling so- 
lution o^^. It has been shown in jsl] that the EdS model 
{yt^i^^q^^ = (1, 1/2) is in addition a saddle point for the 
dynamics of a universe model described by the averaged 
equations. Even a small initial backreaction - which is 
always present due to the observed emergence of struc- 
ture - will drive the expansion away from it and will 
result in accelerated expansion if the deviation goes in 
the corresponding sector, which is the one on the lower 
left-hand side. Besides the EdS model, further special 
solutions are the line r = 1/3 (r has been introduced 
in Eq. ([27|) ). corresponding to models with Friedman- 
nian kinematics but rescaled cosmological parameters, 
and r = representing models without backreaction. In 
this case curvature reduces to a constant curvature 
behavior and they are therefore scale-dependent Fried- 
mannian models. The line r = —1/3 comprises models 
for which backreaction acts like a scale-dependent cos- 
mological constant. The introduction of a genuine cos- 
mological constant would simply shift the diagram down 
by its value if it is positive. 

Figure [71 shows the form of our solutions in this space. 
The dashed curve is the one for the f-regions, the dotted 



one is for M.. Both begin at the EdS model in the center 
as they are initially matter dominated. Because of their 
-limit they lie on the line corresponding to r = —1/5. 
Their opposite signs that are responsible for the vanish- 
ing initial backreaction on V make them evolve from the 
center in different directions. Both lines get shallower 
and £ approaches the line with Friedmannian kinematics 
(r = 1/3) until structure formation sets in and the grow- 
ing variance of expansion rates drives it to the line of a 
scale-dependent cosmological constant (r = —1/3). The 
deceleration parameter is for M not the best quantity as 
it is not guaranteed that it does not diverge for djr = 0. 



V. POWER SERIES OF THE BACKREACTION 
TERMS 

A. Calculation 

In Section |IV] we have seen that the backreaction and 
curvature terms show an initial a^^-behavior. In the fol- 
lowing we will explore where this comes from and rederive 
this behavior in our present context of the partitioning 
into two subregions. We first present the result in the 
form of a proposition: 

Proposition For every evolving spacetime that pos- 
sesses a flat Einstein-de Sitter limit at the beginning of 
its evolution, the backreaction and curvature dependence 
on the volume scale factor ajr on arbitrary subregions 
T stemming from a metric-compatible partition of com- 
pact domains in the hypersurfaces of constant time V 
and that are evolving differently from one another, may 
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be expressed as a Laurent series beginning with aj} for 
backreaction and with for curvature. 

To prove this proposition we will have to find the ex- 
pressions for Qjr (ajr) and {IZ) j, (aj^). They can be de- 
rived from dll) and ([5]) which lead to 



Qm 



ttM (t) 



3 a%, (t) 



(38) 
(39) 



where k = [\ajy^VL^H'^^Y . To calculate Qjr (ajr) and 
iJZ)-p{aj^) we therefore have to derive aM{ciM) ^ind 
{clm)' We achieve this by the assumption of a flat 
matter dominated initial universe model that leads to 



av (t) kt^ ; av {t) = - 



2 ^3/2 2 k^ 



(40) 

To connect to ax) we use the relation aM = clv^m 
and expand with a small parameter a around a = 
which gives 



(41) 



as 



where we have defined \mi and A^v^^ 

'= A^(0)/(3Aa4J and A^^^ •= 
{-2X% (0) + 3A^ (0) Xm.) I (18A^J. A prime stands 
for a derivative with respect to the argument and (0) 
is denoted as A^. := Xm (0). Differentiating (|4T]) with 
respect to time gives (d^^) and aM (d2^,dx»). Using 
(f4Q|) we find clm {av) and d^ (ai^), which leads to the 
finally necessary clm {(^m) and o^m {^m) by an inversion 
of (jil]). This provides 

Qa^ = ?AX,^^'(3A^^+7A^,)— (42) 



(7^) 



40, S ,3, a 



M 



(43) 



and therefore proves the proposition if we take into ac- 
count that the expansion (|^T]) may be extended to arbi- 
trary order. The third- and fourth-order terms of Qm 
are 

-h^ (2A^^ - 3Xm,Xm^ - 9Xms) o? (44) 

2k^ (3A^^ — OAxs^A^^ — ^Xm^^Xmi 

+A^2 + 11Aa^4) X^^a^aM , (45) 

and all higher-order terms may be calculated in a 
straightforward way, but have increasingly complicated 



coefficients. By Equation ((2]), {aj:) and {lZ)j,{aj^) 
are linked and the terms a J of (7^)jr {clj^) are the 
same as those for Qjr (ajr) simply multiplied by = 
— (n + 6) (n + 2)~^. For the second-order term this may 
be seen in (|42]) and (|43]) where = —5. For the third- 
and fourth-order terms above one arrives at {IZ) j, (ajr) 
when multiplying the third-order term with = —3 
and the fourth-order one with = —7/3. 

These results show why we had to include the condition 
of a different evolution of the subregions. If they expand 
in the same manner, the ratio Xm stays constant which 
means that Xm^ = and therefore Qjr = = 0- An- 

other way to phrase it is, that there is no backreaction 
in a completely homogeneous universe. But if there are 
structures developing at any scale, we will have a non- 
vanishing backreaction term and it will go as in a 
matter dominated era. 



B. Quantitative conclusions 

Two remarks are in order here. First we want to em- 
phasize that the result of Equations (|^ and (|^ c oin- 
cides with what has been found by Li and Schwarz(26l,[28| 
in second-order cosmic perturbation theory. They also 
found that on a region P on a flat matter dominated 
background the backreaction term is of second order and 
scales as d'j-^ . The third-order term is a (cosmological) 
constant. 

Second, we want to explore what the result means 
quantitatively for the initial backreaction. For simplicity 
let us first take the fit to the TV-body simulation data of 
Section HVl For the functional form (l36l) we obtain 



«A4 = ^M^av I 1 



1 



2 2 



(46) 



and we may derive the expansion coefficients to be A^^ = 
and Xm^ — —1/3. The first observation is that Xmx — 
causes the term of the averaged scalar curvature 
(7^)jr to be zero. In view of the definition of A^^ below 
Eq. (pij) it becomes clear that this is the case whenever 
^M (^) ~ ^- This means that, if the difference between 
the subregions at the beginning is not very important so 
that they evolve initially in a similar way, the term that 
evolves as a constant curvature term in a FLRW model 
with will be vanishing, provided the background does 
not have a constant curvature term by itself. 

If we demand that the initial matter densities are sim- 
ilar on all the subregions = we find to leading 
order 
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(47) 



which gives the initial value of the backreaction com- 
pared to the initial matter density. This relation is in- 
teresting since it encodes the relation between the evo- 
lution of structure in the Universe and the initial value 
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redshiftz scale R in /z ^Mpc 

Figure 8: Left: Plot of the evolution of the volume percentage of the overdense regions Xm (z) for scales of = 
1, 5 , 10, 15, 20, 30, 50h~^Mpc (from bottom to top). The data points derived by the mesh method of Appendix IB II from 
|60l | are fitted with the functional ansatz of (|48|) . Right: double-logarithmic plot of the scale-dependence of the magnitude 
of the a^-term in the Laurent series of Qm (|42p at the initial time ti (the scale R has been divided by two to match the 
widely used top-hat scales). To illustrate the trend of the data points we added a line with ^-scaling. The ^-axis shows the 
quantity ~ ^q^% while the x-axis shows the top-hat grid scale in h~^Mpc. 



of the backreaction term. In (|36|) a := aj^ character- 
izes the rate at which the Al-regions decouple from the 
common evolution of the subregions. For a larger a, Xm 
decreases more rapidly. This faster evolution requires a 
bigger initial amount of backreaction as it encodes the in- 
homogeneities that force the system to deviate from the 
initial near to homogeneous state. The larger the inho- 
mogeneities are, the faster this deviation takes place. 

To give a quantitative estimate we will however use a 
slightly different fitting function than (|36|) that is adopted 
to the shape of Xm {clv) on a larger range of scales. We 
analyze to this end the A/'-body data provided by (60| 
with a mesh of varying grid size. The result are the data 
points of the left graph of Fig. [51 

They are shown together with a fit using the two pa- 
rameter model 

Y 1 + yaaT)) J 

where 7 (a, j3) is determined by the requirement that 
{(^v) is flat for ax) = I- The form of this function 
is motivated in Appendix IB II 

An interesting outcome from this figure is that it con- 
firms the idea that Xm niay be a good parameter charac- 
terizing the formation of structure. This is because the 
different scales evolve differently as expected from the 
hierarchical formation of structure. On small scales, e.g. 
for R = lh~^Mpc, Xm has already dropped to half of its 
initial value at a redshift of approximately z = 6. For 
scales of R = 10h~^Mpc this happens only dit z ^ 0.7. 
Therefore, we recover at least qualitatively the fact that 
structures start forming at small scales and only then 
assemble to bigger ones. 

Using an expansion of the form (|^6]) and the general 



form of dm) 

n-^' 1 

7j:[- = -7(3A^.+7A^.)«', (49) 

we finally derive the initial abundances of the aj^ back- 
reaction term that are shown on the right-hand side 
of Fig. [HI We have evaluated grid lengths from 1 — 
50h~^Mpc. The line of Fig. [8] indicating an ^-behavior 
of the scale dependence shows that in our evaluation 
of the A^-body structure formation there seems to be 
a discrepancy with results from j23|, where the authors 
found that the values for the backreaction term as de- 
rived from the power spectrum should scale as R~^ with 
the length scale R. In their case, however, they used a 
pure Harrison-Zel'dovich spectrum scaling as to de- 
rive this result. For a more realistic CDM power spec- 
trum, the small-scale behavior is rather k~'^ which in 
turn means that the spatial dependence changes from 
to R~^. The region that we could evaluate with 
scales up to 25h~^Mpc (top hat), is in this small-scale 
regime. It would be interesting to evaluate a bigger sim- 
ulation to go beyond this scale and to see the change in 
the scaling behavior of the backreaction term. 

Quantitatively, the percentage of initial backreaction 
shown in Fig. [8] ranges from 10~^ — 10~^ depending on 
the scale (see Fig. [8]). These specific values show that 
backreaction is indeed perturbatively small in the quasi- 
homogeneous epoch of the Universe. On large scales it 
seems insufficient to change the overall behavior as indi- 
cated by the value of 2 x 10~^ on the scale of the hori- 
zon, found for example in the perturbative calculation 
of [41l|. On the scale of structures, however, this per- 
turbative contribution grows with respect to the matter 
content due to its -scaling and is sufficient to lead to 
the diverse structures we see in today's Universe. 
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VI. CONCLUDING REMARKS AND OUTLOOK 

Let us discuss what we have seen in this paper. Sec- 
tion [TTl showed that a consistent spht of the dynamical 
equations governing the averaged universe model is pos- 
sible. This split also allows, in Appendix [Aj to shed 
light on the property of the volume scale factor to show 
accelerated expansion, even in a case where there was 
no acceleration in the evolution of its components. The 
split of the equations enabled us to construct an averaged 
universe model without having to know the initial mag- 
nitude of backreaction. Instead this latter follows from 
the strength of structure formation that we put in in the 
form of today's volume fraction of the initially over- 
dense regions A4. This model implies that the volume 
scale factor aj:, evolves quite similarly to the scale factor 
in a flat Friedmann model with about 70% cosmologi- 
cal constant. An obvious shortcoming of this model is 
a mismatch with the actual time-evolution of the best- 
fit volume fraction with that predicted by the A^-body 
data. Surprisingly, this mismatch disappears neatly, if 
we reinterpret the fundamental Dark Matter fraction in 
the matter parameter by the backreaction effect, leaving 
only the baryon content as fundamental. We do not want, 
however, to emphasize this result, since the Dark Mat- 
ter issue has to be investigated much beyond our simple 
model. 

Having found that there may be more to the evolu- 
tion of Qjr than just extrapolating the perturbative a'^^- 
behavior, we investigated in Sec.HVlhow to connect struc- 
ture formation and accelerated expansion within a strong 
backreaction scenario. By globally imposing the particu- 
lar scaling solution that mimics a cosmological constant 
and by assuming the structure formation history of a 
Newtonian A-body simulation, the initial a^^-limit on 
the subdomains is obtained, but this also provided a non- 
perturbative extension of the Qjr-scaling to later times. 
The emergence of the a^^-scaling was then closer inves- 
tigated in Sec. El where it became clear that it is generic 
for any matter dominated universe model that starts out 
close to homogeneity. This also confirmed that our choice 
to invoke the -behavior on the evolution on Ai and 
£ instead of V was the right one. We could instead have 
identified our P-region with the one that Li and Schwarz 
looked at, but Section IVl showed that our P-region cor- 
responds rather to their background. Finally, we showed 
that if the Universe may be described by the average 
equations and we require it to have the structure we see 
today, this means that the backreaction component has 
to be of the order of 10~^ — 10~^ (as compared to the 
matter density and depending on the scale one is con- 
sidering) in the initial near to homogeneous phase that 
may be treated by perturbation theory. In the later non- 
perturbative stages of the evolution of the Universe, this 
tiny initial fraction is sufficient to give rise to structures 
and inhomogeneities that we see in the recent epoch. 

A striking result of this work is the fact that both com- 
plementary scenarios lead to qualitatively similar evolu- 



tion laws for the backreaction terms on the largest scales, 
while the models only differ in the details, e.g. in the 
concrete form of the structure formation histories and in 
the behavior of the variables on Al-regions. Since the 
assumptions underlying these two scenarios are quite or- 
thogonal, we are confident that the -scaling behavior 
at early stages of backreaction evolution on the subdo- 
mains (imposed in the first scenario and derived from 
the latter) should be close to the actual physical evolu- 
tion that has to be confirmed by future investigations of 
perturbation theory on an evolving background and by 
nonperturbative models that contain exact solutions as 
limiting cases. 

One out of a number of problems that remain is to 
give a clearcut interpretation of the quantities calculated 
in terms of observables. As we argued it is difficult to 
link the initially underdense ^-regions directly to today's 
voids as one might wish to. To bypass this ambiguity one 
could imagine to extend the analysis to three regions Ai , 
A, and £ where one could choose which overdensity they 
possess. It would be natural to consider as -regions the 
ones with over densities of typical clusters, the f -regions 
with those of voids and put all the astrophysically un- 
spectacular rest into the ^-regions. The advantage of 
probably being able to give a more refined meaning to 
for example the l^-parameters is, however, counterbal- 
anced by the inconvenience to have to introduce another 
parameter in the scaling model or even a free function 
Qa {(^a) in the general case. First results along the lines 
of the scaling of Section IIIII and the fit of Section IIVI 
have been obtained in (50| , but further analysis would be 
needed to find out whether the identification is possible 
and if it is therefore worthwhile to go this way. 

A further refinement of the presented model would be 
to give up the preservation of the identity of the ini- 
tially characterized subregions. Whether this is neces- 
sary could be checked through a detailed analysis of N- 
body simulations. As our model is effective one could 
think of implementing reaction rates between over- and 
underdense regions in the spirit of nonequilibrium transi- 
tions in chemical reactions. A mass-action law could be 
devised that governs an equilibrium state between ele- 
mentary subregions in which the reaction rates are equal 
but nonvanishing, and the transition laws could be de- 
termined as done in nonequilibrium chemical reactions 
(details of such an approach in the cosmological context 
may be found in (65|). 

We shall, however, first investigate other models that 
are based on relativist ic Lagrangian perturbations, gen- 
eralizing the Newtonian nonperturbative model investi- 
gated in m, [131 • This systematic attempt - that is cur- 
rently in preparation - will provide the first generic rel- 
ativistic evolution model for structure formation, includ- 
ing the spacetime metric for implementing observables as 
measured along the light cone. 
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Appendix A: Accelerated expansion in averaged 
models 



We want to use this appendix to elaborate on the pos- 
sibility of accelerated expansion with local deceleration. 
This property of averaged models often leads to confu- 
sion so it may be worthwhile to give a simple example, in 
which one can understand intuitively its origin, and then 
recap some arguments of the literature for the general 
case. One popular paper that gave rise to the confusion 
is, for example, [6g|- Apart from the acceleration issue 
discussed in this section, it also postulates the negligibil- 
ity of backreaction effects due to the applicability of the 
post-Newtonian metric. The reasons of why this argu- 
ment is false, are by n ow given in various papers in the 
literature, e.g. [H, [g^, [Hil . 
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Figure 9: Toy model of a two component model to illustrate 
accelerated volume expansion. On the left the initial config- 
uration with a matterless void 8 and a massive shell M. The 
void in the middle expands in the sequel due to its negative 
curvature. The shell does not expand any more as it contains 
only a virialized matter configuration. This leads to the final 
state depicted on the right where the void volume is dominat- 
ing and the matter regions have been stretched and decreased 
in thickness. In between, the expansion of the total radius a-v 
experiences an acceleration as is argued in the text. 



1. A toy model 

Let us consider the setup of Fig. [9l where we have 
identified two distinct regions JVl and 8. 

The spherical ^-region in the middle is assumed to be 
devoid of matter and to expand only due to its negative 
curvature. The shell however is assumed to be matter 
dominated and virialized such that the net expansion rate 
Hm is zero. This means that the volume of the M region 
is constant, i.e. clm = 0, and that clm = 0- Let us 
consider the case when dg is also vanishing which leads 
with (dni) to 

— =2XM{l-XM)Hj . (Al) 
ax) 

As long as the inner region keeps expanding we will also 
have acceleration of the global expansion, even if none of 
the two regions themselves accelerate. In our case this 
can be understood as a volume effect. As a^: = 0, d^: is 
constant and therefore the growth in radius within a time 
interval At is Aag and also constant. That this does not 
result in a constant expansion rate for the total radius ax) 
is due to the condition that the volume of Ai is constant. 
With an increasing inner radius, the shell must become 
thinner to satisfy this condition. In a more intuitive way 
one may think of the growth of the inner radius as if 
it grew into the shell so that we had to subtract a shell 
with thickness Aag from the inner border of the A^-shell 
and to add the missing volume at the outer border of the 
whole configuration. As the radius of the shell to add 
there is bigger than at the inner border, we need not to 
add Aas but a shell of a smaller thickness to keep \M\ 
constant. This means that the radius ax> of the total 
configuration does not grow by Aag but by Aax) < Aag . 
With increasing size of the expanding inner sphere and 
the corresponding reduction of the thickness of the M- 



shell, as shown in Fig. on the right, the difference of 
the thickness of the shell taken away at the inner border 
of A4 and added at the outer border decreases. Aax) be- 
comes more and more Aag, i.e. it increases. Therefore, 
a-z) experiences accelerated expansion as the growth Aax) 
in a fixed interval At is growing. From this explanation 
it becomes also clear that accelerated expansion can only 
be a temporary effect. In Eq. (|Aip this is reflected by the 
fact that with the growth of the f-region the parameter 
Xm is going to zero and with it dx)- This explanation also 
shows that the effect depends on the number of dimen- 
sions in which one is evaluating it. In a one-dimensional 
setup, it will not appear at all as the two scale factors 
just add linearly. It then grows with every dimension one 
adds what may be seen in the growth of the second term 
in (p!9|) . Adding generic volumes of the same shape but 
evolving differently, one has geometric effects on the evo- 
lution of the global scale factor (e.g. combine 4 squares 
to a bigger square but, as soon as you evolve the indi- 
vidual squares differently, the overall shape will not be a 
square) . 

If we generalize the setup a little to include a decel- 
erated expansion of the ^-region, Eq. ()Aip acquires an 
extra dg term. If it is smaller than the Hg term, one has 
acceleration of the whole configuration despite decelera- 
tion of the inner constituent. 

Note that the acceleration described above also shows 
up in a configuration where the inner part is contracting, 
i.e. Hg < 0. Then, the difference between the shell that 
has to be taken away at the exterior of the outer M- 
region and the one that has to be added on the inner side 
of the A^-shell will be growing. First, they are similar, 
but finally, a reduction of the diameter of the f-region 
by Aag will result in a smaller reduction of the overall 
diameter ax>^ i-e. Aax < Aag. This is to be interpreted 
as acceleration of the P-region in the sense of shrinking 
deceleration. 

To go one step further one could imagine to have a toy 
model of the Universe build up out of balls of the first 
type. Each of them expands in an accelerated way. If 
the accelerated phase sets in during a narrow interval of 
time for all of the balls, one would find that the effect 
adds up and can cause global acceleration. The setup 
in this extended case is a bit similar to the Swiss cheese 
models, but without the restriction of the embedding in 
an overall FLRW evolution. 



2. More general cases 

The emergence of accelerated volume expansion de- 
spite local deceleration has been already widely discussed 
in the literature. Essentially it is due to the fact that the 
local contributions are correlated in an average To 
give a more intuitive explanation, Rasanen, for example, 
argued in [l6| that the physical reason for this volume 
acceleration is that the volume fraction of the faster ex- 
panding regions rises. In our partitioning approach this 
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is reflected in Eq. ([IHl), which may be recast in the form 
aj^ aM as 9 

(A2) 

It underhnes the statement that there can be accelerated 
expansion if the change of the volume fraction of the 
faster expanding regions Xm is sufficient. It also shows 
that the possibility of acceleration holds beyond the toy 
model case. If structure formation is rapid enough so 
that Xm is able to counterbalance I and dg/ag^ it 
will drive acceleration. 

From the averaged equations themselves one can de- 
rive the condition for an acceleration of the volume scale 
factor. For strong backreaction (See [67| for the distinc- 
tion between weak and strong backreaction) Equation 
dll tells us that Qt) must be positive and bigger than 
47rG {g)j) to lead to acceleration. In view of the deflni- 

tion dH]) this means that the term {(O'^)^ — (0)'^^ has to 

be sufficiently bigger than the shear term (cr^)^^, i.e. the 
Universe should be dominated by expansion fluctuations 
rather than by fluctuations of the averaged rate of shear. 



Appendix B: Evaluation of the A^— body simulation 

Here we want to describe how the results on the evolu- 
tion of the Xm parameter in Figs. [4land[5]were obtained. 
The data that we analyzed were obtained from a simula- 
tion of the Virgo Supercomputing Consortium [60]. They 
trace the structure formation in a cube of a sidelength of 
479/i~^Mpc. The underlying cosmological model was a 
ACDM model with I^a = 0.7 and Ho = 70kms-^Mpc-^ 
The simulation contained 512^ ^ 134Mio particles. 



1. A simple mesh method 

To obtain a rough approximation of the overdensity 
fleld, we separated the data cube with a mesh with a 
flxed grid size and determined the number of data points 
in the cells. Then we added up the volume of the most 
dense cells until the sum of the points contained in the 
cells added reached half the number of total points of the 
simulation volume. This is because we decided to flx the 
characterization of the volumes in the initial, near to ho- 
mogeneous state. Under the assumption of a Gaussian 
distribution of the density fleld with only small density 
fluctuations, approximately one half of the mass in the 
Universe will be in over dense M- resp. underdense E- 
regions, if their typical size was nearly the same at that 
epoch. The volume obtained by the addition described 
was suspected to be Vm and used to determine Xm by 
dividing by the total box volume. Of course this identifl- 
cation is a source of error because it may happen that in 
a region that was underdense in the beginning structure 
formation leads to a density peak that belongs to the 



densest half of the matter distribution. But as the con- 
siderations in this paper are mostly concerned with the 
influence of structure formation on the global expansion 
history, we expect the error to be tolerable. 

Apart from this problem the determination is well- 
deflned even if the A^-regions scale differently than the 
global aACDM-scaling of the simulation. This is because 
we are only interested in volume fractions and assume our 
background region V to scale as aACDM in Section HV} 
even if we assume different physical reasons for this be- 
havior. Then 

— — ^M^box , initial a ACBM _ (Bl) 

Vd Nj^Vbox, initial a j^Q-^y^ 

where Nm and Nx) are the number of boxes of the over- 
dense half and the total number of boxes respectively, 
and the simulation volume just drops out. 

To derive a more reliable estimate for Xm than the one 
obtained by the procedure described above, one would 
have to consider Lagrangian domains as required by the 
averaging formalism. This would imply to have a clear 
definition of the volume associated to the particles in the 
initial state and to follow exactly this volume through- 
out the evolution. As in the real world, at least at small 
scales vorticity is not absent, it would be complicated to 
put this into practice, especially in a relativist ic frame- 
work. Other complications include particle crossing and 
domain merging. We suppose, however, that the outcome 
would not be completely different, since averaging strate- 
gies to find extensive quantities over smoothed-out sin- 
gularities are implementable and would likely not change 
the partitioning. It is clear that there is a huge potential 
to improve on our rough evaluation. 

The above procedure was carried out for different grid 
sizes to get an impression of the dependence of on the 
grid length. The result is shown on the left-hand side of 
Fig. [HI We used values of l-50h~^Mpc. The plot makes 
clear that there is in fact a strong dependence of Xm on 
the grid size. This is not surprising as structure forma- 
tion in the Universe proceeds in a hierarchical way. The 
small-scale structure forms first and then starts to com- 
bine to larger structures. This is refiected in the steeper 
decrease of for smaller grid lengths. 

To motivate the choices of grid resolution in Sec. IIIII 
we consider that a scale of lh~^Mpc, as used for the 
lowest curve of Fig. [TOl may already be too small as the 
simulation is used for large-scale structure. The highest 
curves above 10h~^Mpc grid length is on the other hand 
already merging typical overdense structures like clusters 
with underdense voids. If one takes into account that a 
typical cluster of galaxies is in the range of 2 — 8Mpc one 
would expect that grid sizes of 2 — 8h~^Mpc are perhaps 
best adopted to the considered problem. Therefore, we 
chose a grid length of 5h~^Mpc for the comparison of 
Fig. SI 

For the results of Section IV Bl we had to find an ade- 
quate functional form for the evolution of Xm {av)- From 
the theoretical point of view it should meet several re- 
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redshift z logarithm of the density contrast +1 



Figure 10: Left: Comparison of the results for a box size of = 1, 5, 10, 15, 20, 30, 50h"^Mpc (from bottom to top) with 
the one obtained by the Voronoi approach (dots on the curve that represents the fit of Fig. [5]). The corresponding curve varies 
more strongly especially at the beginning which may be interpreted as faster structure formation. Right: distribution of the 
density contrast that has been assigned to the simulation points by the Voronoi method. The lowest curve shows the initial 
distribution of z = 5. The other curves show the same distribution at redshifts of 3, 2, 1, 0.5 and (from bottom to top). 
The values on the x— axis represent log (^ + 1), the y— axis the number of regions with the corresponding density contrast in a 
logarithmic scaling. 



quirements. First of all, we expect to start at an ini- 
tial value of about 0.5 on all scales, as the density fluctu- 
ations of the Early Universe are assumed to be Gaussian. 
The second condition should be that it does not deviate 
rapidly from this value at the beginning. So Xm should 
have a flat tangent at = 1. Third, we expect that 
tends to on all scales as matter clusters strongly and 
the voids keep expanding. For a standard ACDM model 
this is even more true as then the cosmological constant 
will accelerate the growth of the empty regions between 
the matter filaments. These three conditions lead to the 
ansatz of (l48l) . 

2. The Voronoi method and structure formation 

As the fixed mesh is not adapted to the real matter dis- 
tribution in the simulation volume, we used in a second 
step a partitioning of the volume into Voronoi regions. 
This approach has been advocated by van de Weygaert 
(HI, [t^ and assigns each simulated point the region that 
is closer to it than to any other point. In condensed mat- 
ter physics a Voronoi cell is better known as the Wigner 
Seitz cell. To calculate the Voronoi cells we used the 
program "qhull" [71], and derived Xm in the same way as 



described above by adding the densest Voronoi cells. The 
resulting points are shown in Fig. [10] on the left together 
with the curves for 0.5, 1 and 5h~^Mpc, as derived by the 
mesh method. The plot shows that the Voronoi results 
are closest to the lh~^Mpc-mesh ones. The shape, how- 
ever, shows a steeper behavior in the beginning, which 
one would interpret as a faster structure formation. The 
Voronoi result was used for the fit of A^ of Fig. [5l 

We finally use the separation into Voronoi cells to show 
the formation of structure in the distribution of the over- 
density field. 

To this end we plot in Fig. [10] on the right the evo- 
lution of the distribution of the density contrasts of the 
Voronoi cells. For our first data at a redshift of z = 5 
the distribution still is strongly peaked around zero. The 
maximum is then moving to under dense regions. In the 
meantime a second maximum emerges for high density 
regions. The mean over density, however, stays zero all 
the time. It is interesting to see that the distribution 
spreads and reaches an extension of 8 orders of magni- 
tude in the density contrast. These extreme differences in 
the local densities and their expected different evolution 
is one of the main motivation for the idea of a structure 
formation effect on the global evolution. 



